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Abstract 

We establish a decoupling result for the P and S waves of linear, isotropic elasticity, in the setting of twice- 
differentiable Lame parameters. Precisely, we show that the P O S components of the wave propagation 
operator are regularizing of order one on L 2 data, by establishing the diagonalization of the elastic 
system modulo a L 2 -bounded operator. Effecting the diagonalization in the setting of twice-differentiable 
coefficients depends upon the symbol of the conjugation operator having a particular structure. 

1 Introduction 

We consider the linear, isotropic elastic wave equation on R™, and Cauchy initial value problem on a time 
interval [—to, to), 

dfw l -(A(x,D)w) l = 0, 

Wi\ t =0 = fi, (1.1) 

d t Wi\ t =o = gt, 
where w — (wi, . . . , w n ) denotes the displacement vector, and 



n 

(A(x,D)w)i =^2d k (n(x)(d k Wi + diw k )) + 6\( X(x)d k w k ) . (1.2) 

k=l 



The elastic wave equation arises in many physical situations, including the propagation of seismic waves 
through the Earth. The identification and separation of P and 5* polarized constituents (so-called phases) 
underlies many applications of seismic data processing and imaging, which traditionally has been based on 
the assumption of smoothness of the medium and asymptotic considerations. As an example, we mention the 
notion of receiver functions, which were introduced and developed by Vinnik (1977) and Langston (1979). 

Throughout this paper, we assume that the Lame parameter functions /i and A belong to C 1 ' 1 {W l ), in 
that the first derivatives V/i and VA are globally bounded, Lipschitz functions on R™. 

The principal symbol <j2{A)(x,£) of A(x,D) takes the form 

-a 2 (A) = fa) (CI 2 I + (Ms) + A ( a; )) £ ® f 

(1.3) 

= (2fi{x) + X(x)) |e| 2 <7 (IIp) + p(x) ICI 2 a (U s ) 



1 



where ITp and lis are the projection operators 

Upw = A _1 V(divi(;) , H s w = (I - U P )w 
corresponding to the Fourier multipliers 



The symbols of lip and lis are thus, for each x, £, projections onto the eigenspaces of A(x, £), respectively 
of dimension 1 and n — 1. This yields that 

A(x. D) Upw = (2fj,[x) + \(x))AU P w + l.o.t. , 

(1.4) 

A(x, D) Tl s w = (J,(x)A U s w + l.o.t. , 

where l.o.t. denotes an operator of first order in w (on some finite range of Sobolev spaces, as discussed in 
the next section.) 

We call Upw and II5W, respectively, the P and 5* modes of the elastic wave w. If the Lame parameters are 
constant, then the wave propagation operator commutes with lip and II5, so that wave evolution preserves 
modes which are of pure P or pure S form. 

In the case of nonconstant \i and A, there will in general be a nonvanishing coupling of the P and S 
modes in the evolution. The main result of this paper is that the P O- S terms in the wave propagator are 
one order smoother than the P — » P and S — > S terms. Precisely, let the solution to (|1.1[> be written in the 
form 

w(t,-)=c{t)f + s(t)g, 

and decompose c(i) = cpp(t) + cps(t) + csp(t) + css{t), where cpp(t) = IIpc(t)IIp, etc. We similarly 
decompose s(t). The operators cpp(i) and css{t) are bounded on a range of Sobolev spaces H s (M. n ), in 
particular for < s < 1, and spp(i), sss(t) map ff s_1 (t n ) to H s (M. n ) for the same range. (All maps are 
strongly continuous in t.) We prove that 

c PS (t) : L 2 (R n ) -> iJ^R") , s PS (t) : ff- 1 (R") -> ff^R") . 

Additionally, we show that cpp(t) agrees with the propagator for App = ilpAUp up to a similarly regular- 
izing term; see Corollarv l3.3l for a precise statement. 

Establishing regularization for the P -O- S components of the propagator requires diagonalizing A(x, D) 
to second order; that is, one order beyond the principal symbol consideration of (|1.4j) . since treating a first 
order derivative of w as a driving force term will produce terms of the same order as w in the solution. 
This diagonalization is carried out through conjugation by I + K, where K is a compound pseudodifferential 
operator with finite regularity symbol of order —1. 

In the case of smooth \x and A, a similar diagonalization was carried out, in a more general setting, to 
all orders in Taylor (1975) using the pseudodifferential calculus. He used this diagonalization to study the 
reflection of singularities for the system of isotropic elasticity, as well as for other systems with constant 
multiplicities. Stolk and de Hoop (2002) have used Taylor's diagonalization in the study of the linearized 
seismic inverse scattering problem. Dencker (1982) analyzed the propagation of polarization in Euclidean 
space for linear systems of principal type with smooth coefficients, which includes isotropic elasticity. Wang 
(1998) applied Dencker's propagation of polarization result to the linearized seismic inverse problem; in 
particular he discusses the issue of polarization in reflections caused by discontinuities of a perturbation of 
the medium. 

In the case of C ' coefficients, the construction of K is more delicate, and carrying through the compo- 
sition calculus to first order depends on the symbol of K having a special form. We believe that this is the 
first article studying the interaction of P and S waves for the elastic system in non-smooth media, and it re- 
mains open under what conditions this can be carried out in the C > setting for more general diagonalizable 
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systems with constant multiplicity eigenvalues. We anticipate that our result will have applications in the 
study of the linearized seismic inverse problem in media of limited smoothness. There has been considerable 
activity in studying unique continuation in the stationary case for non-smooth Lame parameters, see Lin 
et. al. (2010) and the references therein for the most recent results. 



Notation. 

All pseudodifferential operators that arise in this paper will have symbols that are finite sums p(x, £) = 
Pm-j(x, 0) where p m -j is homogeneous in £ of degree m — j. We use a m (jp) to denote the principal 
symbol p m (x, £) of p, where m is the order of p. 

We will use L 2 (R n ) and H s (M. n ) to denote the L 2 and L 2 -based Sobolev spaces, in the case of complex 
valued functions as well as n-dimensional vector valued functions, such as the solution w. 

Throughout the paper, R will denote an operator that is bounded on L 2 (R ra ); the form of R may change 
in each occurrence. 



2 The Conjugation Operator 

In this section we construct the conjugation operator K, which diagonalizes A(x, D) up to a L 2 -bounded 
operator. We decompose A = App + Aps + A$p + Ass, where 

Ap P = n P A(x, D) n P , A PS = n P A{x,D)Ils, 

Asp = H s A(x, D) n P , A ss = n s A(x, D) W s . 

Note that Asp = A PS , while both Ass and App are self-adjoint. Furthermore Asp, and hence Aps, is of 
order 1, in that 

A PS : H s+1 (R n ) -> H s (W n ) for - 2 < s < 1 . (2.1) 

To see this, we write 

A PS w = A' 1 V(div A n s w) = A' 1 V( [div, A] TL s w) , (2.2) 

where we use div lis = 0- Each component of the composition then takes the form A -1 V(d(b dUgw)), where 
b denotes a first-order derivative of fi or A, which being Lipschitz acts as a multiplier on H s for — 1 < s < 1. 
This establishes (|2.1j) for such s. If instead one writes 

A PS w = —Up A (A - V div) A -1 w 

and commutes one derivative from the right past A, one obtains (|2.1[) for the range —2 < s < 0. 
Our main conjugation result is the following. 

Lemma 2.1 There exists an operator K of order —1, in the sense that 

K : H s (R n ) -4- H s+1 {K n ) for - 1 < s < 1 , (2.3) 

such that 

U P (I + K)A(x, D)w = App {I + K)w + R P w , (2.4) 

U S (I + K)A{x, D)w = A ss (I + K)w + R s w , (2.5) 

where Rp and Rs are bounded maps on L 2 (R"). 
Additionally, 

Kpp = K SS = , K SP = -K PS , (2.6) 
where Kpp = lip Klip, Ksp — Us KTip, etc. 
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To establish (|2.4p ~ (|2.5l) . we expand A and K in terms of their P and S components. Terms such as UpKsp 
and AppKsp vanish. By (12.11) and (|2.3|) . the terms KpsAsp and KspAps are bounded on L 2 (R n ), hence 
are absorbed into R. Using (|2.6|) . the conditions (|2.4[) - (|2.5|) thus reduce to the following conditions 

KpsAss - AppKps + A PS = R, (2.7) 



KspApp - A ss K S p + A SP = R, (2.8) 

meaning precisely that the left hand sides are bounded maps on L 2 (R"). Condition flU) is equivalent to 
(|2.7|) by taking adjoints, so we will establish (|2.7j) . 

Consider first the case of smooth ^ and A. Since Kps is of order —1, and App and Ass of order 2, the 
equality (|2.7[) reduces to saying that the top (i.e. first) order terms in the symbol expansion for the left hand 
side cancel out. By (|1.3p . this implies that 

r, \ -(Tl(Aps) , . 

(MO) + A(x))|£| 2 

We will also reduce ()2.7j) to a principal symbol calculation for the case of // and A of regularity C 1 ' . Care 
needs to be taken, however, since by (|2.2j) the symbol of Aps, hence Kps, is only Lipschitz in x, whereas 
App and Ass are of second order. That we can reduce matters to a symbol calculation depends critically 
on the fact that Kps has a particular structure, which in turn depends on the special form of Aps- 

We consider the operator Aps more closely. Since divng = 0, by (|1.2[) we can write 



(A PS w)i = ^2 A 1 d i d j d k (n (dk(ttsw)j + dj(n s w) k ) 



= 2 ]T A^didj ( (dkn) dj(Usw)k) 

j,k=l 

n n 
k=l j,k=l 

In going from the first to the second line we used that J^k 9k(^sw)k = 0, together with symmetry of the 
expression in j, k. The second term in the third row is a bounded operator on L 2 ; hence we may write 

A PS w = 2 V( (V/i) • U s w) + R . 

Motivated by (|2.9p . we now define 

K PS w = 2<S>{D)A- 1 V((fi + A)- 1 (V A i) -U s w) , 
and set Ksp = —K* PS so that (|2.6[) holds. Here, $(^) is a smooth cutoff to |£| > M, 



\(\>2M, 
\Z\<M, 



with M depending on the C 1,1 norms of [i and A. Condition ()2.3[) holds since (A + /i) 1 Vfi is a multiplier 
on fP for — 1 < s < 1. Furthermore, by taking M large, we will have 

\\Kw\\ H . < %\\w\\ H . , -1<s<1. (2.10) 

Finally, we note that since Kpsw is the gradient of a function, it follows that 

K PS = npifps = Kp S Tls ■ 
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To evaluate the composition AppKps, we observe that we may write 

n 

App = (2(i + A) AIIp + R k dk , (2.11) 

fe=i 

where 

R k : H s (R n ) -> iT a (R n ) , -1<s<1. 

This follows from the fact that the commutator [Updk, fj] is bounded on H s (M. n ) over this range, which is a 
consequence of the Calderon commutator theorem; see Calderon (1965). 
Since d k ° Kpg is bounded on L 2 , it follows that 

AppKpsw = 2 (2(i + A)$(D)V( (/x + A) _1 (Vit) • n s w) + Rw 

— 2 <1>(_D)V( (2/i + A) (/x + A) _1 (V^) • U s w) + Rw . 

To evaluate the composition KpsAss, we first observe that we can write 

Assw = iIsA(/xn5w) + A\W , 

where A\ : L 2 — s> if -1 . To see this, we note that div lis = 0, so the A terms in A do not contribute to Ass- 
For the same reason we can write 

Iis^dk((Ldi(Usw)k) = Us y^(dfc^) dt(U s w)k 

k fc 

which maps L 2 — > H . Finally, we write 

UsY.dki^dk^sw),) =U s A((iU s w) l -U s J2(dk^dk(U s w) l , 

k k 

and observe that the last term maps L 2 — > H^ 1 . 

By (|2.3j) . the composition Kp$Ai is bounded on L 2 , hence we are reduced to showing that 

$(L»)A- 1 V( ((i + A)- 1 (V/i) • (A - Vdiv)(^n s u;)) = $(£>)V( (fi + A)- 1 (V/i) • (fiU s w)) + Rw , (2.12) 

where we write IlgA = A — V div, and understand that A acts on each component of a vector, i.e. as a 
scalar operator. 

The relation (|2.12|) is derived as a consequence of the identities 

A(V- W) = V- AW + 2 div (W- W) - W- (AV) (2.13) 

V- V div W = div (Vdiv W) - div (Wdhi V) + W-VdivV (2.14) 

where in the term div (W ■ VV) the div pairs with V. 

Setting V = (fi + A) _1 V/x, and W = (iHsiv, then (|2.12p is equivalent to 

$(L>)A" 1 v( (V- (A - Vdiv)VF - A(V- W) \ = Rw . 

By (|2.13I) - (I2.14I) we can write the left hand side as 

$(£>)A -1 v(W- (A - V &xv)V - div (V&xvW) + div(WdivV) -2div(W- VV)) . 

The last two terms, 

$(£>) A -1 V div( Wdiv V - 2 W- W) 
5 



are L 2 bounded operators applied to W, hence of the form Rw, since V is Lipschitz. Since divW^ = 
(V/x) • Hsw, the same holds for $(Z?)A -1 V div(Fdiv W). The remaining term is where we use the special 
structure, since it involves two derivatives of V. The key observation is that 

(A- Vdiv)V^ = (A- Vdiv)((// + A)- 1 V/i) € L°°(M n ) , 

since AV/i = VdivV/i, and hence the expression in fact involves only second derivatives of fi and A. Conse- 
quently, 

w^-^(D)A- 1 v(w-(A-Vdiv)v) : L 2 (R") -> i/^M") . 

3 The Parametrix and Decoupling of Modes 

We will realize the parametrix construction for (jl.ll) using parametrices for the purely polarized equations, 
following the procedure introduced in Smith (1998). 

Theorem 3.1 There exist strongly continuous one parameter families of operators 

c P {t) : H s {W n ) -> H s (R n ) , s P (t) : iI s_1 (R n ) ZP(R n ) , 

d t c P (t) : H s (R n ) -> iF-^R") , 9 t s P (t) : H a ~ 1 (R n ) -4- ff s - x (R n ) , 
/or < s < 2 and £ € [— tojto]) suc/i that the solution to the Cauchy problem 

(<9 2 - App)v P = G P , v P \t=o = fp, d t v P \ t = a = g P , 
where Tlsfp — Tlsgp — lis Gp = 7 and 

/p£ff s (l"), 3P 6i? s - 1 (r), GpeL^Mcfol.F- 1 ^")), 

is given by 

Vp(t, • ) = c P (t)fp + s P (t)gp + / s P (t - s) Gp(s, ■) ds . 

Jo 

The same result holds if P is exchanged with S, and both results hold for all to > 0. 

Proof. Following pages 816-818 of Smith (1998) we produce families of operators c(t) and s(i), respectively 
bounded and regularizing of degree 1 on H s (R n ) for all s, with c(0) = I, d t c(0) = 0, s(0) = 0, <9ts(0) = 7, 
such that for < s < 3, 

(5 2 -(2 M + A)A)c(t) =T (t) : H s (R n ) -> fl < - 1 (R n ) , 

(3.1) 

(a 2 - (2/i + A)A)s(t) = r x (t) : ff s - x (R n ) -> Jf s - x (R n ) . 
We set cp(t) = IIpc(t)rip, Sp(t) — IIps(i)rip, and observe that, for < s < 2, 

(a 2 - App)c P (t) =T Pfi {t) : H s (R n ) -> ff-^R") , 

(3.2) 

(a 2 - App)sp(t) = r P ,i(t) : F"- 1 (R n ) -> iF^R") . 
To see this, we use the decomposition (|2.1ip to see that 

Appcp(t) - (2/i + A)Acp(t) : iJ s (R") -> Jf^^R") , < s < 2 . 
As in the discussion following (|2.11l) . the commutator 

[Hp, (2/x + A)A] : H s (W l ) -> iP-^K") , < s < 2, 
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and (|3.2p follows from (I3.1[) . The exact evolution operators cp and sp are obtained from cp and sp following 
pages 819-820 of Smith (1998), using a convergent iteration. For example, 

cp(t)=c P {t)+ [ s P (t- s)T{s)ds, 
Jo 

where T(t) is the solution to 

T(t)+ f T Pjl {t-s)T(s)ds = -T Pfi (t). 
Jo 

We observe here that, since cp(t), sp(t) and Tpo(i) are orthogonal to Us on both the right and left, it 
follows that so is cp(t). 

The bounds on dtCp(t), etc., follow similarly, using the fact that dtc(t) and dts(t) respectively map 
H s (R n ) -s- H a - l (R n ) and ip-^R") -> i/^^R") for all s. □ 

Suppose that w is the solution to the linear elastic system (jl.ip . We pose 

v = (7 + i^)ui , wp = ITpw , = lis?; ■ 

By (|2.4[) we can then write 

<9 t 2 wp = lip (J + iT)A(a:, D)w 
= A PP (I + K)w + Rw 
= Appvp + R(I + K)^v , 

where at the last stage we use (|2.10j) to invert (I + K) on the range of Sobolev spaces we will consider for 
w. Similar considerations apply to vs, hence we have the system 

(dj - App)v P = R P v , vp\ t=0 = fp, d t vp\ t=0 = g P , 

(3.3) 

(df - A ss )v s = Rsv , v s \ t=0 = fs , d t v s \t=o = gs , 

where Rp and R$ are bounded on L 2 (M"). Here, we have set fp = Hp(I + K)f and gp — Hp(I + K)g, and 
similarly for f s , gs- 

The system (|3.3I) is equivalent to the integral system, 

v P (t) = c P (t)fp +s P (t)gp + / s P (t- s)R P v(s, -)ds, 

Jo 

(3.4) 

vs{t) = c s (t)f s + s s (t)g s + / s s (t - s)R s v(s, -)ds. 



We now restrict attention to data / e H s (R n ), g e H"- l (R n ), with < s < 1. Then 

cp(t)fp + s P (t)gp e c(Ho,t ],ff s (r)), 

and the same holds with P replaced by S. For < s < 1, then sp(t — s)Rp is a bounded map on H s (M. n ), 
and similarly with P replaced by S, so the Volterra type system (|3.4[) can be solved by iteration to yield a 
solution 

v = vp + vseC([-t ,t Q lH s (W n )). 

The solution takes the form 

v = c P (t)fp + s P (t)g P + c s (t)f s + s s (t)gs + r(t)(f,g) , 
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where r(i) : L 2 (R") © ff- x (R n ) i2" x (R n ) and t r(i) : L 2 (R") © TJ-^R") -> L 2 (R"), with both r(i) and 
d t r(t) strongly continuous in t. 

We next write f = (I + K)f, g = (I + K)g, and w = (I + K)- X v. By (l2~3l) . 

(I + if) -1 - J : fP(R™) -¥ H s+1 (R n ) for -1<s<1. 

Consequently, 

(/ + X)- 1 c P (i)(7 + /f) - c P (t) : L 2 (R") -> ff x (R n ) . 

By similar results for the other terms, we have the following corollary (where r(i) is of a different form than 
above). 

Theorem 3.2 Let w(t,x) solve the Cauchy problem (fTTTj) . luitfc dato / £ L 2 (R n ), ,g G F _1 (R"). Let 
(fp,9p) = n p(/,3), and (fs,9s) = H-s(f,g)- Then w may be written as 

w(t, • ) = cp(t)f P + B P (t)g P + c s (t)fs + s s (t)g s + r(t)(/, <?) , 

w/iere r(t) is bounded from L 2 (R") © i2"- x (R n ) to and d t r(t) is bounded from L 2 (R n ) © ff- 1 (M n ) 

to L 2 (R"). 5oi/i r(i) and 9jr(t) are strongly continuous in t. 

In particular, if the initial data for tti is a pure P mode, in that /<? = 55 = 0, then 

n sW = n S T(t)(f, 9 ) e c([-t ,*o],ff 1 (R n )) , 

provided that /p e L 2 (R") and gp € iJ _1 (R"). This can be viewed as saying that the P S coupling is 
regularizing of order 1, and the same holds for pure S modes. 

This decoupling can also be expressed in terms of mapping properties of the components for the solution 
operator to (| 1 . 1 [) . An immediate result of Theorem 13.21 is the following. 

Corollary 3.3 Let the solution to (jTTTJ) be written in the form 

w(t,-)=c{t)f + s(t)g, 

and decompose c = cpp + cps + c$p + css> s = s pp + s ps + s sp + s ss, where cpp = rtpcllp, etc. Then 
with cp , cs , sp , sg as in Theorem \3.1[ we have 

cpp(t) - c P (t) : L 2 (R n ) -> i? x (R n ) , 
cp<?(» : L 2 (R n ) ->if x (R n ), 

spp(t) - s P (t) : H-\R n ) -> i? 1 (R n ) , 
sp S (t) :ZT- 1 (R n )^-ff 1 (R n ). 
results hold with P and S interchanged, and are strongly continuous in t. 

Furthermore, in each case the time derivative of the operator is a bounded map from the indicated domain 
into L 2 (R n ). For example, 

d tCpP {t) - d t c P (t) : L 2 (R") -> L 2 (R n ) , 
with strongly continuous dependence on t. 

The corollary follows from Theorem 13.21 since all operators listed arise as compositions of r(t) with the 
projections lip and lis- 

We remark that if the stronger regularity assumption is made, fi, A € C^' 1 with k > 2, then the range 
of Sobolev indexes s on which r(t) is regularizing of order 1 increases, and Corollary 13.31 thus applies to a 
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range of H s (R n ). Alternatively, by adding successive terms to K, the system (I3.3[) for v = (I + K)w can 
be decoupled to higher order, that is, with Rp and Rs of negative order; see Taylor (1975) for the case of 
smooth /i, A where Rp and Rs are smoothing operators. This requires adding lower order terms to App 
and Ass- The order of gain in regularity in Corollary 13.31 cannot be increased, however, since it is limited 
by the term of order —1 in K. 
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